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Abstract 

A simple scaling theory for the sintering of fractal aerogels is presented. The 
densification at small scales is described by an increase of the lower cut-off 
length a accompanied by a decrease of the upper cut-off length £, in order to 
conserve the total mass of the system. Scaling laws are derived which predict 
how a, £ and the specific pore surface area £ should depend on the density p. 
Following the general ideas of the theory, numerical simulations of sintering 
are proposed starting from computer simulations of aerogel structure based on 
a diffusion-limited cluster-cluster aggregation gelling process. The numerical 
results for a, £ and £ as a function of p are discussed according to the initial 
aerogel density. The scaling theory is only fully recovered in the limit of 
very low density where the original values of a and £ are well separated. 
These numerical results are compared with experiments on partially densified 
aerogels. 

PACS numbers: 61.43.Hv, 64.60.Ak, 81.20.Ev. 
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I. INTRODUCTION 



Sintering of silica aerogels is a process in which the original material is heated at a tem- 
perature smaller than the melting temperature of the corresponding silica crystal. It results 
in a strengthening of the internal structure and a gradual elimination of the pores accom- 
panied by a general densification and shrinkage of the whole sample. Such process allows 
the design of intermediate materials, hereafter called PDA (partially densified aerogels) of 
increasing density up to full dense silica glass. Several theories have been introduced to ex- 
plain the shrinkage due to sintering ]T|-0. It is now quite well established that the internal 
mechanism for sintering of aerogels is due to local viscous flows of matter. For the applica- 
tion to glasses and aerogels, the most convincing theoretical approaches are due to Scherer 
@||. Here, we present a scaling approach, which can be considered as a generalization of 
Scherer's type of calculations to fractal matter. Since a short account of the theory has 
already been published elsewhere [[J, we present here a slightly different (simpler) deriva- 
tion of the scaling laws (part 2A). Based on these theoretical ideas, we present two recently 
developed numerical procedures able to simulate the sintering of realistic aerogel structures 
(part 2B). Then, after presenting the numerical results (part 3), we compare them with 
experiments (part 4) and we conclude (part 5). 



II. CONSTRAINTS ON THEORY 



A. Scaling theory 

It is now well established (see for example the recent interpretation of low angle neutron 
scattering experiments 0]) that silica aerogels are made of connected fractal aggregates. 
Fractal scaling occurs in a range of lengths extending from a lower cut-off a, which is 
the mean diameter of the silica particles constituting the aggregates, up to an upper cut- 
off £ which is the mean diameter of the aggregates, or, equivalently, the mean center-to- 
center distance between two neighboring connected aggregates. Since it is admitted that, 



2 



in aerogels, sintering results from local viscous flow of matter, we will assume that only 
the shortest lengths are concerned and therefore only the lower cut-off a is affected, without 
changing neither the fractal dimension of the aggregates nor their mutual arrangement. After 
a first stage where the individual shapes of silica particles are only slightly modified, while 
bridges are built between them, one reaches a "scaling" regime where individual particles 
can no longer be distinguished. In this regime, we will assume that there still exist a lower 
cut-off a, which can then be defined as the mean thickness of the aggregates arms. Since 
the sintering process tends to lower the internal surface pore area by transfering matter 
from large curvature to small curvature regions and therefore reinforcing the thinner arms, 
there results a gradual increase of a. But, correlatively, to insure the mass conservation, 
there should be a general shrinkage of the material, characterized by a decrease of £. Figure 
1 gives a two-dimensional sketch of how an aggregate looks like at two different stages of 
sintering. 

Since we assume that the overall inter-aggregate structure is conserved, the shrinkage at 
the scale of an aggregate is the same as at the macroscopic scale. Therefore, if p is the mean 
density of the material, the following "trivial" scaling should hold: 



Then, using the results of the fractal geometry for simple scale-invariant structures ||, the 
minimum number N of balls of diameter a necessary to cover the total mass of a fractal 
aggregate of diameter £ (the requirement being that each point of silica matter should be 
inside at least one ball), is given by: 



Very often, this formula appears in the literature with only one parameter (a or £) because 
the other one is assumed to stay constant. Given a, the relation N ~ £ D corresponds to 
the "mass-size" relation. Given £, the relation N ~ a~ D corresponds to the Haussdorf- 
Besicovitch definition of the fractal dimension (for simple self-similar fractals all the defini- 
tions of the fractal dimension lead to the same value |§). But here, we need to keep the 



f ~ P 3 



(1) 




(2) 



a 
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two parameters together since they both vary. Using this formula, the bulk density can be 
expressed as a function of a and £: 

Na 3 ,a. 



p~ ^ ~ qr D (3) 

Then, combining with equation (1), we immediately get: 

a ~ fPW-m (4) 

Even if the present derivation of equation (4) is more straightforward than the one already 
given ||, it implicitely uses the same simple symetries (fractal scaling and mass conserva- 
tion). 

As soon as a and £ are known, any other structural property, which is related to these 
parameters, can be calculated. Several examples have been previously provided but here 
we will focus on the specific pore surface area E which is the area of the whole silica-air 
interface counted per unit of silica mass. Following the above ball-covering reasoning, the 
total interface for one aggregate, is of order Na? while its mass is of order Na 3 , therefore 
one has: 

£ „ ^! „ I „ p -D/^-D) (5) 

Na 3 a w 

A key point of our reasoning is that a remains a well defined cut-off. In other words the 
constitutive aggregates can be modelized by mass fractals downs to length a. Therefore to 
satisfy (1) and (3) simultaneously, a should increase if £ decreases and reciprocally. Such 
reasoning fails if the interface might be considered as a surface fractal; then one could imagine 
to smooth the surface without changing £. This may happen for other kind of materials and, 
here, for the last stages of sintering where the scaling theory doesnot work. 

Before pursuing, it might interesting to compare our scaling results with previous 
Scherer's approach of sintering ||||. In its first calculation ||, Scherer was considering 
a regular cubic array whose bonds were made of cylinders of length £ and diameter a. Ne- 
glecting surface deformations near the cylinders connections and simulating the sintering as 



an increase of a/£, he was able to calculate analytically the density as a function of a/£, at 
least for a < £: 

^ = 3vr(f) 2 - 872(f) 3 (6) 
Ps 4 4 

where p^ is the density of silica. This formula can be usefully compared with our formula 
(3). Since Scherer is considering the peculiar case of D = 1 fractals his first term is the same 
as ours, but, in his simple geometry, he is able to perform the full analytical calculation and 
he gets a second term which, in our language, appears to be a correction to scaling since it 
is negligeable compared to the first term when £/a is large. 



B. Numerical simulations 



Since it has been shown that the diffusion-limited cluster-cluster aggregation model p.lC 



provides a quite realistic modelization of an aerogel structure ||?|,|TT|j, in this section we 
intend to use such a model to develop numerical simulations of aerogel sintering. This is 
still an approximate treatment of sintering since we will use some naive coarse-graining 



procedures, such as that used in real-space renormalization group methods jn|, to describe 
the smoothing at the shortest scale, but, since we are now considering a realistic structure, 
we hope to get some informations on corrections to scaling. 

We have developed two numerical methods. The first one works on a cubic lattice and 
considers cubic particles of edge length ao (volume Vq — <Zq, mass m = psafy- The second 
one works off-lattice and considers spherical particles of diameter ao (volume vq = |Oq, mass 
m = Ipsa 3 ,)- In both cases the particles are inside a cubic box of edge length La and their 
number N Q is such that the volume fraction cq is set to a desired value cq = N vo/ L 3 al- 
Therefore, except L which should be chosen as large as possible, the only parameter of the 
model is c which is directly related to the initial aerogel density p — c oPs- 

To build the original (non-sintered) aerogel structure, the particles are first randomly 
disposed (i. e. put on randomly chosen sites, avoiding multiple occupancy, on lattice, or 
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sequentially centered at random points, avoiding overlaps, off-lattice) in the box. Then, 
these particles are allowed to undergo a brownian diffusive motion and they irreversibly 
stick when come on contact. Aggregates of particles are also able to diffuse together with 
the individual particles and to stick to particles or to other aggregates. In this diffusive 
motion, the diffusion constant of the aggregates is considered to vary as the inverse of their 
radius of gyration and periodic boundary conditions are assumed at the box edges PHTT 



When the concentration cq is sufficiently large (larger than a threshold value c g which tends 
to zero for infinite box size f7j), the final structure is a gelling network which extends from 
edge to edge in the box and which can be described as a loose random packing of connected 
fractal aggregates, of fractal dimension D ~ 1.8, whose mean size £ decreases as c increases. 

The modelization of the sintering process differs if one considers the on-lattice or the 
off-lattice version of the model. 

1. On lattice 

Here we have used a discrete blocking method (based on "box counting" ideas [§]) which 
proceeds step by step. It implies that the value of L is set to some power of 2, L = 2 n . At 
step p of sintering, the original cubic lattice is replaced by a "super-lattice" of parameter 
length 2 p ao, each super-cell of the new lattice containing 2 3p cells of the original lattice. 
Then a super-cell is considered as occupied if it contains at least one particle of the original 
structure. The number N p of occupied super-cells is computed. Then the sintered structure 
is obtained by considering the structure made of the N p occupied super-cells and by applying 
an adequate length contraction (3 P in order to conserve the total mass. The volumic fraction 
c p of the super-structure being: 

Np^l_ N p 



the contraction factor is given by: 



= 2^)* (8) 

Co JV 



and the actual value of the cut-offs a p and £ p at step p are given by: 



,0,0 ,N ps 



i 



f, = # (10) 

As expected when combining formulae (8) and (10) one recovers that the "trivial" scaling law 
(1) on £ is automatically satisfied, but the one on a deserves to be tested by the numerical 
calculations. On figure 2 (a) we provide a typical example where a section of the box is 
shown at three different steps of sintering. 

To determine the specific pore surface area, we numerically determine the number S p 
of square interfaces separating occupied supercells nearest neighbor to a non-occupied one 
(taking care of the periodic boundary conditions) and we calculate E by: 

_ 1 S p a 2 p _ 1 S p 
p s N p a 6 p a p s N s N s 

In practice, for a given Co and at each step p, the quantities N p and S p , calculated numerically, 
have been averaged over several realizations of the initial configuration. 

Note that we could have removed the restriction to supercells of edge lengths 2 p a , 
working instead with edge lengths £ao, £ being any integer. Imagining such an extension 
of the method, let us call if the value of £ above which all the supercells become occupied 
and define £o as being £/ao- With this definition £ is not, rigourously speaking, the mean 
aggregate diameter but, rather, it is such that £o — do represents the edge length of the 
largest cubic hole in the structure. Using then the fact that, at the end of the sintering 
process, one has £/ = a/, one gets: 

a f = £ c 1 (12) 

In practice we have written a separate code to evaluate £o by determining if. Therefore 
using (12), one can get a numerical estimate of af. 
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2. Off lattice 



Here we have used a dressing method which, compared to the preceeding method, has 
the advantage to be continuous. We first consider a s-dependent dressed structure, where s 
is a continuous variable, in which each initial sphere of diameter a is replaced by a sphere 
of the same center but of larger diameter a^s), given by: 

a>d(s) = a (l + s) (13) 

Then the total volume V d (s) located inside the overlapping spheres is numerically calculated 
(avoiding multiple counting of overlaps). The volumic fraction of the dressed structure is 
given by: 

Here also, to insure mass conservation, the sintered structure is obtained from the dressed 
structure after applying an adequate length contraction (3(s) given by: 



Then the actual diameter a(s) of the spheres of the sintered structure and the correlation 
length £(s) are given by: 

a(s) = m = ao m (16) 

^ - W) (17) 

Here again the scaling law on ^(s) is automatically verified while a(s) needs a numerical 
calculation. The value for ^ is estimated from a/, the limiting value of a(s) when c(s) 
tends to one, by inverting formula (12). Here £ — a represents the diameter of the largest 
spherical hole. On figure 2(b) we provide a typical example where a section of the box is 
shown at three different steps of sintering. 
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To determine the specific pore surface area S(s), we observe that the total surface Sa(s) 
of the dressed structure is simply related to the derivative of the function Vd(s): 

= 2 ,m = i<^w 

aa^ a as 

Then, dividing by the total mass and correcting by the adequate scaling factor, we get the 
following expression S(s): 

w >t 2 1 dV d a( \ 2 dlogcjs) 
aops V d (s) as a p s as 

In practice, for a given Co, we have numerically determined the whole curve c(s) versus s as 

well as its derivative and we have calculated, a(s) and S(s) using formulae (16) and (19). 

Here also, for each c value, the results have been averaged over several realizations of the 

intial configuration. 



III. RESULTS 

In both the on-lattice and off-lattice cases, we have made calculations with L = 64, we 
have considered initial concentrations c = 0.01,0.02,0.05 and 0.1 and the results for c, 
a and £ have been averaged over 20 independent realizations of the initial configuration. 
Figure 3 gives the results for a/a and £/a as a function of c (log-log plot) for the different 
c values. Cases (a) and (b) correspond to on-lattice (symbols) and off-lattice (continuous 
lines) results, respectively. The slope indicated is the theoretical slope _D/3(3 — D) — 0.5 
expected from the scaling theory with D — 1.8. In principle, the scaling theory should be 
recovered in the vanishing concentration limit c — > where a and £ are well separated. 
Even if the results for a versus c are more and more linear, with an apparent slope close 
to 0.5, when c decreases, in both cases we observe, down to c = 0.02, some corrections to 
scaling both at low c and large c. 

The low-c corrections to scaling are different in the two simulations. When performing a 
low-s expansion in the off-lattice formulae, it can be shown that a in this case should behave 
quadratically when c tends to Co: 
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(^(off-^) 2 (20) 
a Cq 

This singular behavior is due to the peculiar non sintered structure made of perfectly tangent 
spheres. In all other situations, with non zero areas between connected particles, as it is in 
the on-lattice case, the low-c behavior should be linear: 



,a — aos c — Co ( , 

\ jon ~ (21) 

ao Co 

Therefore, in the off-lattice case, there is an initial regime during which finite areas are grown 
between connected particles. During this regime, for the same increase of p, a increases less 
than in the on-lattice case. After this regime, the evolution of a becomes similar to the one 
of the on-lattice case, but for a larger initial concentration. This is shown in figure 4 (a) 
where we have reported the off-lattice results for Co = 0.01 (dashed curve), which, except 
at low concentrations, are very close to the on-lattice results for Co = 0.02. Note that, 
apart from low-c corrections to scaling, the lattice structure built with spherical particles 
of diameter ao, instead of cubes, which has a concentration smaller (c = (7r/6)co), behaves 
quite similarly during sintering. This is reasonable, since, starting from the same spheres 
it is known that the on-lattice and off-lattice diffusion limited cluster-cluster aggregation 
processes lead to the same large distance correlations [|T^ . 



At large c (c close to 1) the correction to scaling comes from the fact that, in this limit, 
it remains only a few holes in the structure and, obviously, one can no longer speak of a 
structure made of connected fractals. The volume fraction of the last hole being 1 — c and 
its radius being proportional to aj — a, where a/ is the limiting value of a, one should have: 



a / -a~(l-c)a (22) 

This behavior is well verified by both on-lattice and off-lattice numerical results. 

Concerning the specific pore surface area, we have calculated the dimension-less quantity 
aoPo^-- It is worth noticing that we get different results already for c = cq. In the off-lattice 
case, we recover the trivial result: 
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(a po£ ) ff = 6 ( 23 ) 

which means that the structure is made of non-overlapping spheres. In the on-lattice case, 
we get a lower value since we do not count the interfaces between neighboring particles. 
If we neglect the contributions of loops, it can be shown that the 6 should be replaced by 
4 — 2/iVo (since there should be Nq — 1 bonds). But, due to the presence of loops, the actual 
number is smaller, close to 3.5 and decreases when Co increases. The numerical results for 
£ as a function of c are given in figures 4 (a) and (b) as a log-log plot of versus c for the 
different c values. In both cases the slope -0.5 expected from the scaling theory is better 
verified than for a. Here also the cq for the off-lattice case is close to the curve for c' = cq/2 
for the on-lattice case (but not as close as in the case of a). When c tends to 1, we recover 
that £ tends to zero, which is an improvement compared to the toy model. Following the 
above reasoning for c close to 1, one should get: 

S ~ — - — ~ (1 - c)~3 (24) 
(if — a 

a behavior quite well observed in our off-lattice results. 

Another result of our numerical calculations is the estimate of £o, which is here related 
to the largest spherical hole in the non sintered aerogel. In table I, we compare our off- 
lattice results with the alternative estimates (here called £q) obtained from the location of 
the minimum of the pair correlation function g(r) 0. It is reasonable that these estimates 
are close and roughly proportional to each other. 



IV. DISCUSSION 

In this section we would like to discuss our numerical results at the light of several 
experimental results on partially densified aerogels (PDA). Let us discuss first some previous 
low angle neutron diffraction experiments |TJ]]. In figure 6, we show the scattering intensity 



I(q) curves for a series of PDA of increasing densities made by sintering a "neutral" aerogel, 
hereafter called N46. These experimental results are in good qualitative agreement with the 
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theory. The slope of the fractal linear regime (here corresponding to D ~ 2.3) does not 
depend on the density, except in the last steps of sintering (where anyway the slope can no 
more be interpreted as a fractal dimension), but its q extension is gradually reduced as p 
increases. The arrows show how the parameter a and £ have been previously estimated, i. e. 
by assuming that a -1 and £ _1 correspond to the lower and uper departures from linearity. 
At the light of the recent interpretation of the I(q) curves for non-sintered aerogels @, we 
now know that such method is approximate and cannot give the right absolute values for 
both a and £. In particular, since £ should be related to the position of the low-g maximum, 
not visible on the figure, its actual value should be considerably larger than the above 
estimate. However, since the estimation has been done with exactly the same procedure for 
all densities, we can think that the relative values a/a and £/£o are relevant. The same 
method has been applied for another neutral sample, N26, of lower initial density but of 
similar fractal dimension and on a basic sample, B46 of lower fractal dimension D ~ 1.8, 
for which it is believed that the DLCA model applies 0. 

Figure 6 (a) gives the results obtained with the basic sample as a log-log plot of both 
a/ao and £/ao versus p/ps and the best fits of a/a^ with both calculations. The continuous 
curve is the result of the off-lattice calculation with c = 0.04 and the circles correspond to 
the on-lattice calculation with Co = 0.08. Since the actual density of the non-sintered sample 
is about 0.2g/cm 3 , the experimental Cq value should be about 0.09, closer to the value taken 
for the on-lattice fit. The fact that the on-lattice model works better is certainly due to the 
presence of finite areas between connected particles already in the non-sintered material. 
Such "pre-sintering" might have occured during the supercooling process used to obtain the 
aerogel from the gel by extracting the solvant. The data for a are quite well fitted except 
for the point of largest density. Anyway, since this point corresponds to a very large density, 
the method of determining a from the I(q) data becomes highly questionable. Note that in 
the last stages of sintering, one cannot interpret I(q) by the fractal theory and one needs 
another numerical calculations of I(q) (such as the scattering by a collection of random 
polydisperse holes of radii a/ — a at averaged mutual distance £) to be able to interpret the 
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scattering curves. Moreover the numerical estimates of £ (except the last point) exhibit a 
nice —1/3 slope providing a strong support to the main hypothesis of our theory. However 
the estimates are smaller (by a factor almost 5) than the theoretical values, as expected if 
one remembers that the theoretical £o value is related to the location of the maximum of 
the I(q) curve. 

Figure 6 (b) gives the experimental results for two neutral samples with different initial 
densities. Here we cannot try to quantitatively compare with our numerical results since 
the fractal dimension is larger. It is known that the formation of neutral aerogels cannot be 
explained by the DLCA model. They have a tenuous and flexible polymeric structure and 
therefore deformations and restructuring cannot be avoided during aggregation. However the 
effective slope of the a versus c curves is larger than in the basic case in perfect aggreement 
with the scaling theory, as already mentionned in ||. Assuming a fractal dimension of 
2.3, the theoretical slope D/3(3 — D) should be slightly larger than 1. Assuming that the 
corrections to scaling enters similarly than in the basic case, the experimental slope is smaller 
as seen in the figure. Moreover, here again, the results for £ exhibit a nice —1/3 slope. 

V. CONCLUSION 

In this contribution we have presented a scaling model for the sintering of aerogels, we 
have etayed this approach by some numerical calculations and we have compared the numer- 
ical results with some experiments. The numerical calculations have shown that corrections 
to scaling are not negligeable, in particular at the last stages of sintering where the fractal 
scaling process is no more valid (the notion of cut-off fails) and where the sintering process 
rather involves a gradual elimination of residual holes. In the first stages of sintering the 
corrections to scaling have a larger role in the c-dependence of the lower cut-off a than in 
that of the specific pore surface area. The effective slope of the numerical log a versus log c 
curve is slightly smaller than the one expected from the scaling theory. These theoretical 
results are in quantitatively and qualitatively good agreement with experimental results on 
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basic and neutral aerogel, respectively. A quantitative agreement in the case of neutral aero- 
gels requires a specific modelization of these materials which have a larger fractal dimension 
than basic aerogels. Nevertheless, the numerical calculations presented here are still using 
crude modelizations of what might be the actual sintering process. Even if the two numer- 
ical methods presented here give almost the same results, in both cases, at each step of 
sintering, the internal surface of the pores is not a minimum area surface as it should be. 
In particular, in the last stages of sintering, the remaining holes are not spherical. A more 
sophisticated numerical approach would be to consider a minimum area surface at each step 
of restructuring, its area decreasing gradually during the sintering process. A Monte Carlo 
modelization of such process is under progress. Moreover, we are presently considering the 
off-lattice numerical approach as a basis to simulate the gas transport properties of partially 
densified aerogels in order to interpret some new experiments on the permeability of gas 
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(Venezuela). 
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TABLES 



c 


a f 


£o 


Co 


0.01 


6.2 


23 


28.8 


0.02 


4.6 


13.5 


16.9 


0.05 


3.2 


6.1 


8.7 


0.077 


2.5 


4.0 


5.9 


0.1 


2.3 


3.0 


5.0 



TABLE I. For each concentration cq considered in the off-lattice simulations, we have 
reported the value of a/, the corresponding value of £ an d the location £ of the minimum 
of the pair correlation function g(r). 
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FIGURES 



FIG. 1. Two dimensional sketch of an individual aggregate at two different stages of sintering. 

FIG. 2. Section of the box for the original configuration and at two stages of sintering, (a) and 
(b) correspond to the on-lattice (with L = 32) and off-lattice (with L = 40) method respectively. 

FIG. 3. Log-log plot of a/ao and £/ao versus c = p/ps for L = 64 and different initial con- 
centrations Co- Cases (a) and (b) correspond to the on-lattice and off-lattice method respec- 
tively. In case (a) open circles, black squares, open diamonds and black triangles corresopond to 
Co = 0.01,0.02,0.05 and 0.1, respectively. The solid lines are the results for £/ao- The dashed 
curve shown in (a) corresponds to the off-lattice results with co = 0.01. 

FIG. 4. Log-log plot of S/So versus c = p/ps for L = 64 and different initial concentrations 
co- Cases (a) and (b) correspond to the on-lattice and off-lattice method respectively. In case (a) 
open circles, black squares, open diamonds and black triangles corresopond to co = 0.01,0.02,0.05 
and 0.1, respectively. The dashed curve shown in (a) corresponds to the off-lattice results with 
c = 0.01. 

FIG. 5. Low angle neutron scattering intensity curves I(q) for a series of neutral PDA with 
increasing density. 

FIG. 6. Experimental dependances of a/ao (black symbols) and £/ao (open symbols) as a 
function of c = p/ps (log-log plots) for a basic sample (case a) and for two neutral samples (case 
b). In case (a) the full line is the theoretical curve obtained with the off-lattice model with co = 0.04 
and the open circles are the results obtained with the on-lattice model with cq = 0.08. 
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fig 3a 




fig 3b 




fig 6b 




